Current sheet bifurcation and collapse in electron magnetohydrodynamics 
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Inertial effects in nonlinear magnetic reconnection are studied within the context of 2D 
electron magnetohydrodynamics (EMHD) with resistive and viscous dissipation. Families of 
nonlinear solutions for relevant current sheet parameters are predicted and confirmed nu- 
merically in all regimes of interest. Electron inertia becomes important for current sheet 
thicknesses 5 below the inertial length d e . In this case, in the absence of electron viscosity, 
the sheet thickness experiences a nonlinear collapse. Viscosity regularizes solutions at small 
scales. Transition from resistive to viscous regimes shows a nontrivial dependence on resis- 
tivity and viscosity, featuring a hysteresis bifurcation. In all accessible regimes, the nonlinear 
reconnection rate is found to be explicitly independent of the electron inertia and dissipation 
coefficients. 
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Magnetic reconnection is a fundamental mechanism for magnetic energy release in both astro- 
physical and laboratory plasmas. It manifests itself as a topological rearrangement of the magnetic 
field lines, followed by a conversion of magnetic energy into particle energy, plasma kinetic energy 
and heat, and is characterized by the presence of localized current sheets. A long-standing problem 
in the theory of reconnection is to identify the relevant microscopic mechanisms that render the 
process efficient, and to predict the transition from slow [as in resistive magnetohydrodynamics 



(MHD)] to fast reconnection 



Pi 



Two-fluid effects enable fast reconnection [2j in MHD. Ions and electrons can decouple in their 



relative motion within some relevant microscopic scale, allowing for enhancec 



the context of the well-known Hall MHD two-fluid model, various numerical 



reconnection rates. In 
[5] and theoretical 



J HE 



efforts have concluded that the transition from slow to fast reconnection occurs when 



di > A/y^S^, where di = cjio v i is the ion inertial length, A is the characteristic length in the 
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Figure 1: Diffusion region geometry. 

plasma outflow direction, and S v ^ 1 is the resistive Lundquist number. 

However, qualitative differences between theory and computations remain. In particular, numer- 
ical evidence of a hysteretic bifurcation, in the transition between resistive and Hall MHD regimes, 



[j], [sl- Des 

fj, E|j an explanation for such strongly nonlinear behavior has remained e^ 

In this Letter, we extend the nonlinear analysis put forth in Refs 
electron inertia. We identify the interplay between electron inertial effects and dissipation as the 
root of the observed hysteretic behavior. For simplicity, we restrict our analysis to the electron 
magnetohydrodynamics (EMHD) model and use it as a paradigm of the more general Hall MHD 
model. In EMHD, the magnetic field is frozen into the electron fluid, while ions are a neutralizing 
background at rest (v, « 0) within the ion inertial scale length dj. Using the stagnation-point 



usive. 

by including finite 



in 



configuration proposed in Refs. 
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111 ], we describe the 2D diffusion region by replacing the full 
EMHD partial differential equations with a low-dimensional dynamical system (few time-dependent 
ODEs) and study its steady-state properties. In such a way, we derive families of nonlinear solutions 
for the diffusion region aspect ratio and the associated reconnection rates. 

Nonlinear Reduced Model. Expressing the EMHD equations in Alfvenic units with di as the 
equilibrium scale length, and retaining electron inertia corrections, gives Q]: 

d t B* + V x (j x B*) = -r?V x (V x B) + rj H V x (V x V 2 B) , (1) 

where B* = B + d^V x (V x B). Here, r\ and tjh are the dimensionless resistivity and electron 
viscosity (or hyper-resistivity), and d e = \J m e /m,i is the electron inertial scale length. 
Writing Eq. ([T]) in component form, assuming d z = 0, gives (with V = i]—r] H V 2 ) 

d t B* - V • (j p B* x - B* p j x ) = -V{d 2 yx B y - dyB x ), (2) 
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d t B* y - V • (i p B* y - B* p j y ) = -V(d 2 yx B x - d 2 B y ), (3) 



8 t B* + d\ (j p • V) V 2 B Z + B p ■ Vj z = VV 2 B Z . (4) 



Here, j p = Vx (B z z) = -v e , B* = (B*, B y), and B* = B z - d 2 e V 2 B z . 

To proceed, we follow Refs. [8|, and consider a rectangular 2D reconnection (diffusion) region 
of dimensions 5 and w (Figure [I]) . We define the upstream and downstream magnetic fields as 
B x = x • B(0,<5/2) and B y = y • B(tt;/2,0), and define the discrete flow stream function B z = 
— z • H(w/2,5/2). Consequently, the inflow and outflow velocities are given by v V) e = —2B z /w 
and v Xie = 2B z /5, respectively. Then, we discretize Eqs. ([2])-(|4|) at (x,y) = (0,5/2), (w/2,0), and 
(w/2,5/2), respectively. Using d x ~ 1/w, d y ~ 1/5 to find B-Vj 2 ~ -(B x /w+B y /5)(B y /w-B x /5), 
and (j p • V)V 2 -Bz ~ [j? — -jjjr] , we obtain a set of equations for B x , B y , and B z (dropping tildes 
and numerical factors of order unity for simplicity): 

a* <5 R * B Z B* \( B v B x\ fK \ 

B * ~ 5 x ~ ~5w~ = m W) {5^-W)' (5) 



w B z B* 

b;--b;+ y 

y w y 
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\ ow w z 



(6) 




B x + B y \/B y _B x \ = / l_ + ±\ B+ 4 B 2 ( J_ _ 1 

w 5 J \w 5 J ' \5 2 w 2 J 5w z \w 2 5 2 

(7) 

where V(6, w) = 7] + rj H (5~ 2 +w~ 2 ), B* = B x + d 2 (B x /5 2 - B y /5w), B* = B y + d 2 (B y /w 2 - B x /5w), 
B* = B z + d 2 (5~ 2 + w~ 2 )B z , and the overdot denotes time derivative. 

Steady-state Solutions and Reconnection rates. Fixed points of Eqs. (El) -(3) provide insight 
into the intrinsic limitations of reconnection rates at nonlinear saturation 



12|. Setting time 



derivatives to zero, and introducing the parameters d e = and £ = — , we obtain from Eqs. J5|) 
and © B y /B x = £ (1 + 2d 2 e )/(l + 2^£ 2 ) , and B Z /^2B X = S' 1 ^ 1 ~ 0/ [l + d 2 e (l + £ 2 )] . Here 
S- 1 = S" 1 + S^it' 2 + 1) is the inverse of the effective Lundquist number, with S v = \f2B x /r\ 
and Sh = V2B x w 2 /t\h- Using these relations in Eq. (jTj) , gives the equation for the diffusion region 
aspect ratio £(S, d e ): 



(8) 




In the massless electron limit (d e = 0), Eq. |(SJ) recovers solutions obtained in Ref. 

The reconnection rate, defined as the electric field in the ignorable direction at the X-point 
(x = y = in Fig. [T]), is given by E z = Vj z \x , where j z \x = (B x /5 — B y /w) is the current density. 
Using the previous results in this expression for the reconnection rate gives: 

R 2 P~ x - P 

E^V^S- 1 ^^ (9) 

w i + 2d? e e 

From Eq. j9]) it is evident that, for given B x and w, large electric fields E z preferentially occur for 
£ 2 <C 1. We consider this limit next. For simplicity, we also assume 2d 2 !; 2 = 2(d e /w) 2 <C 1, which 
is true for small enough d e . Then, since (1 + 2d 2 ) /{I + d 2 ) « 0(1) for any d e , Eqs. ((Hj) and ((9]) 
simplify to become 

e-s- 1 — U^, (10) 
1 + d 2 ' 

E z ^V2^(l + df). (11) 
Fiscows regime (t)h > 0, r] = 0). Rewriting £ = d e /(d e w) in Eq. (TT0|) gives 

I + i"Ujv5^ = '"" (12) 

which implies 5/d e ~ (^) 1/3 for 4 < 0(1) (magnetized regime), and <5/<i e ~ 77^ for d e > 0(1) 
(inertial regime). These scalings have been numerically validated and will be discussed later in this 
Letter. In particular ^n the magnetized regime, 5 ~ (j]uw j \[2B X ) X I^ > d e . In the inertial regime, 
as discussed in Ref. [12|], the plasma is demagnetized within the inertial scale and the bulk current 
thickness is determined by d e , so that j z \x ~ 2B x /5 « 2B%jd e , where B% = ~k- B(0, d e /2) is the 



magnetic field upstream of the inertial region. Then, 5 « \l tihw / {s/2B%d e ) < d e describes the 
radius of curvature of the current sheet at x = 0, which sets the reconnection rate. Employing these 
expressions for 5 in Eq. JTTJ gives for the reconnection rate: 

E ^^USE. (13) 

It! 

where B XiTnax = max^, is the magnetic field at the upstream boundary of the induced current 
i.. Note that the —tion rate in the vi S con S reg i me is not an expHcit Action of eiectron 
viscosity [8] or inertia [13] and is therefore potentially fast. This result implies that electron physics 
is enabling fast reconnection, while, as already suggested in Ref. [ijj], ion inertia can eventually 
limit it (in fact, for arbitrary d { and in the Hall MHD regime, E? al1 « E z di Q, Q, Q). Unlike 
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the massless case d e = 0, electron inertia limits the electron outflow velocities at the inertial scale 
length d e by v x « B z /5 ~ B x /5 < B x ^ max /d e = V^.e, the electron Alfven speed, as expected. 
Resistive regime (rj > 0, rjn = 0)- Rewriting £ = d e /(d e w) in Eq. (fTUI) gives 



Equation ([14]) features a saddle-node bifurcation with a threshold in the parameter rf , such that 
steady-state solutions for c£ e (or 5) exist only for rj* > 2. In the magnetized regime [d e < 0(1)], we 
find a single solution 5 = rjw / (y/2B x ) > d e [8|, and the reconnection rate is given by Eq. ([TO . with 
d e — ► 0. As in the viscous regime, the electron outflow velocity is limited by the electron Alfven 
speed. In the inertial regime [d e > 0(1)], we find 5 « \/lB x d\J {r\w) which, after substituting 
B x = B x -^, results in d e = i]w/(V2B x ) for any 5 < d e . Thus, the quantity 5 is not determined and 
can reach arbitrarily small values below d e . This is a consequence of the fact that small resistivities 
cannot set a dissipative length scale when inertia is important. Indeed, if we introduce ^(x,y,t) 
such that Bp = z x V*, then Eq. p} gives [J f t (* - d 2 e j z ) = rjj z , with d/dt = d t + v e • V, and 

r V 4 

j z = V 2 ^. When * < d 2 e j z , i.e. 6 < d e , we find J z (a;, y, i) ~ 0. This is a hyperbolic equation 

for j Zl which cannot set a dissipative scale, thus it cannot prevent the collapse of the current sheet 
thickness to zero below d e . This result implies that, in the resistive regime, the reconnecting system 
will experience a loss of equilibrium when the parameter rj* becomes sufficiently small, resulting in 
a transition to another state. The nature of this new state critically depends on whether viscosity 
is present or not. 

Hysteresis bifurcation. Equations (fT2l) and (fT4l) are valid in the asymptotic limits rj = and 
rj* H = 0, respectively. The general steady-state solution for the current sheet thickness, for finite 77, 
r]H, and d e , is obtained from Eq. (fT0j) as 5 3 — rj*5 2 + j 2 5 — (irj* H = 0, with 5 = 5/d e . Here, we have 
introduced the empirical coefficients 7 and (3 to take into account multiplicative numerical factors of 



. This equation is known as the universal unfolding 
la ]. It can be shown that the equilibrium manifold 



0(1) neglected in the derivation of Eqs. ©-((I 
of the pitchfork bifurcation of codimension 2 
features hysteresis for rf H < (7/\/3) 3 //3. 

Numerical Validation. We employ the magnetic island coalescence instability to validate pre- 
dictions of the model. The ideal-MHD-unstable equilibrium is given by the magnetic flux function 
^(x,y,t) = —A log [cosh (j) + e cos (%)] [a], where A = 1/2-7T is the equilibrium characteristic 
length scale, and e = 0.2 is the island width. Results are obtained by performing a series of non- 
linear 2D simulations Q] varying rj, rjn, and d e . Values for 5, w, B x are measured at the instant of 
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Figure 2: S/d e at the time of maximum reconnection rate as a function of rj* H = w/ (\/2B x d^) from nonlinear 
simulations. The transition to the inertial regime S/d e < 1 happens at S/d e ~ 0(1) and rf H ~ 0(1). 



maximum reconnection rate, when the process saturates non-linearly, and a current sheet is already 
formed between the two coalescing islands (at y = along the x-direction). The downstream length 
w is evaluated at the point of maximum outflow, B x is measured upstream at (0,5/2), and the 
current sheet thickness 5 = 2^J2 log 2 is found as the full width at half maximum, where y* is 

(l2| . 



defined from d?, j z 

In the viscosity-dominated regime, the scalings from Eq. (fl~2l must hold, and this is what we find 
numerically. In Fig. [21 we show 5/d e plotted against the normalized viscosity rf H = r]Hw/(V2B x d^,) 
for r] = 10 -5 , 7.63 x 1CT 7 < rj H < 7.63 x 10 -6 , and 5 x 1(T 3 < d e < 2.25 x 10" 2 . Both scalings 
S/d e ~ (iJh) 1 ^ 3 f° r S/de > C(l) and 5/d e ~ ry^ for 5/d e < 0(1) are identified, and the transition 
occurs at 6/d e ~ 0(1) and 77^ ~ 0(1), as expected. Numerically, we find 0.55 < E^ um < 0.62, 
which agrees with the prediction of Eq. (fTTI) within a factor of two. 

In the resistive regime, Eq. ([14]) predicts the absence of a steady-state solution for values of 
resistivity such that "if < 0(1). In this case, numerical simulations indicate that the current 
density develops an arbitrarily thin sub-d e nonlinear scale, as shown in Fig. 02 A similar behavior, 
conjectured by Wesson [lfj], was first understood in the framework of nonlinear collisionless tearing 



modes 



17|]. When the threshold condition rf > 0(1) holds, two resistive steady-state nonlinear 



solutions for current layers are found for a certain range of rj* (and one otherwise). The black 
dots in Fig. |4] are resistive results from nonlinear simulations with rjjj = 0, r? = 5 x 10~ 2 , and 
10~ 2 < d e < 2.25 x 10~ 2 . The dashed line is the solution of Eq. fUJ rewritten for 5 = 5/d e , 
25 = rj* ± y (rj*) 2 — 47 2 , where numerically we find 7 = 1.65. In this case, we observe that 
the minimum value 5 m i n ~ d e for the numerically obtained current thickness is always such that 
8 > d e , as explained before. When 5 > d e holds, two steady-state solutions for 5 are possible for 
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Figure 3: Example of nonlinear collapse of the current sheet in the resistive regime. Here 77 = 5 x 10 
di = 1, r] H = 0, and d e = 2.25 x 1(T 2 w 1836- 1 / 2 . 

6/d e 




Figure 4: Current sheet thickness 5 = S/d e at the time of maximum reconnection rate as a function of 
rj* = wr//(V2B x d e ). The dots are for T]h = 0, 77 = 5 x 10~ 2 , and 10~ 2 < d e < 2.25 x 10~ 2 . The crosses are 
for d e = 2.25 x 1CT 2 , and 10" 2 < 77 < 7 x 10" 2 , with r) H = 2.29 x The squares are for d e = 2.5 x 10" 2 , 

and 10~ 2 < 7] < 5 x 10~ 2 , with t]h = 10 -7 . The dashed, solid, and dotted lines are solutions of Eq. lfT0|) . 

S 3 - ifS 2 + 7 2 <5 = Prix, with 7 = 1.65, (3 = 9, and r)* H = rf 'i] H / '(rye? 2 ) . 

27 < 77* < 1 + 7 2 , and the quantity 5/d e is not single- valued in rf (see Fig. H]). 

The presence of electron viscosity regularizes the current density for 5 <d e , and r]* < 2, allowing 
for a nonlinear steady-state current profile with finite thickness. As explained earlier, the transition 
between resistive and viscous regimes is nontrivial. It depends strongly on 77^, and exhibits hysteresis 
for rf H w 3.6 x 10~ 2 < / [3 ~ 0.11 (squares in Fig. [4]), and the lack thereof for rf H w 

2.4 x 10 _1 > 0.11 (crosses in Fig. |4}. We note that the numerical solution seems to be able to map 
all branches of the S-curve. This is likely due to the fact that the island coalescence problem is 
highly dynamic, and the system survives a very short time at the point of maximum reconnection 
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rate. A careful study of the stability properties of the bifurcated equilibrium manifold is left for 
future work. 

In conclusion, we have extended recent steady-state nonlinear reconnection theory to include 
the effect of electron inertia and to study its interplay with dissipation parameters. In the absence of 
electron viscosity, for sufficiently small resistivities, we have confirmed earlier observations of current 



sheet collapse 



16 



171 ] and provided, for the first time, a nonlinear threshold for such behavior. 
Electron viscosity regularizes the current layer at small scales and allows the system to achieve a 
nonlinear steady-state in both inertial (<5 < d e ) and magnetized (5 > d e ) regimes. The transition 
from resistive to viscous regimes shows a nontrivial dependence on resistivity and viscosity. For 
sufficiently small viscosities and for a range of resistivities, three different states are available for 
5/d e (see Fig. [4}. Thus, we conclude that electron physics is responsible for earlier numerical 
evidence of hysteresis jsl, [J. We note that this fact may have been obscured by unrealistically small 
di/d e ratios employed in previous simulations. Finally, in all accessible regimes, the maximum 
reconnection rate is formally independent of electron inertia and both dissipation coefficients. 
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